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Abstract 

This paper contains an analysis of rank-k solutions in terms of Riemann invariants, 
obtained from interrelations between two concepts, that of the symmetry reduction 
method and of the generalized method of characteristics for first order quasilinear hy- 
perbolic systems of PDEs in many dimensions. A variant of the conditional symmetry 
method for obtaining this type of solutions is proposed. A Lie module of vector fields, 
which are symmetries of an overdetermined system defined by the initial system of equa- 
tions and certain first order differential constraints, is constructed. It is shown that this 
overdetermined system admits rank-k solutions expressible in terms of Riemann invari- 
ants. Finally, examples of applications of the proposed approach to the fluid dynamics 
equations in (k+1) dimensions are discussed in detail. Several new soliton-like solutions 
(among them kinks, bumps and multiple wave solutions) have been obtained. 



Resume 

Dans cet article, nous presentons une analyse des solutions de rang k exprimees en 
termes d'invariants de Riemann obtenues a partir des relations entre les methodes de 
reduction par symetries et des caracteristiques generalises pour les systemes hyper- 
boliques et quasilineaires du premier ordre. Une variante de la methode des symetries 
conditionnelles pour obtenir ce genre de solutions est presentee. Un module de Lie de 
champs de vecteurs, representant des symmetries du systeme surdetermine constitue 
du systeme d'equations initial et certaines contraintes differentielles du premier ordre, 
est construit. II est demontre que ce systeme surdetermine admet des solutions de 
rang k exprimees en termes d'invariants de Riemann. Finalement, a titre d'exemple, 
une discussion detaillee de Papplication de l'approche proposee aux equations de la 
dynamique des fluides en (k+1) dimensions est donnee. Plusieurs nouvelles solutions 
de type algebriques, rationnelles et solitoniques (bumps, kinks et ondes multiples) ont 
ete obtenues. 
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1 Introduction 



The general properties of nonlinear systems of PDEs in many dimensions and techniques for 
obtaining their exact solutions remain essential subjects of investigation in modern mathe- 
matics. In the case of hyperbolic systems, the oldest, and still useful, approach to this subject 
has been the method of characteristics which originated from the work of G. Monge [32] • In 
its modern form it is described e.g. in [TT],p6j,[30j,[40j,[41J,[44j. More recently, the devel- 
opment of group theoretical methods, based on the work of S. Lie [29], has led to progress 
in this area, delivering new efficient techniques. However, these two theoretical approaches 
have remained disconnected and have provided, in most cases, different sets of solutions. The 
symmetry reduction methods (SRM) certainly have a broader range of application, while the 
generalized method of characteristics (GMC), though limited to nonelliptic systems, has been 
more successful in producing wave and multiple wave solutions. Thus, the mutual relation 
between these two methods is a matter of interest and we have undertaken this subject with 
the view of combining the strengths of both of them. 

The approach to constructing rank-k solutions presented in this paper evolved from our 
earlier work [22] , [21] , [25] , aimed at obtaining Riemann k- waves by means of the conditional 
symmetry method (CSM). The main idea here has been to select the supplementary differen- 
tial constraints (DCs), employed by this method, in such a way that they ensure the existence 
of solutions expressible in terms of Riemann invariants. Interestingly, as we show later, these 
constraints prove to be less restrictive than the conditions required by the GMC. As a result, 
we obtain larger classes of solutions than the class of Riemann k-waves obtainable through 
the GMC. 

The organization of this paper is as follows. Section 2 gives a brief account of the general- 
ized method of characteristics for first order quasilinear hyperbolic systems of PDEs in many 
independent and dependent variables. A geometric formulation of the Riemann k-wave prob- 
lem is presented there. In Section 3 we reformulate this problem or rather, more generally, 
a problem of rank-k solutions expressible in terms of Riemann invariants, in the language of 
group theoretical approach. The necessary and sufficient conditions for obtaining this type 
of solutions are determined after an analysis of their group properties. A new version of 
the conditional symmetry method for construction of these solutions is proposed. Sections 
4 to 7 present an application of the developed approach to the equations describing an ideal 
nonstationary isentropic compressible fluid. We find rank-1 as well as rank-2 and rank-3 
solutions admitted by the system, among them several new types of soliton-like solutions 
including kinks, bumps and snoidal waves. In Section 8, we construct rank-k solutions for 
the isentropic flow with sound velocity dependent on time only. We show that the general 
integral of a Cauchy problem for this system depends on k arbitrary functions of k variables. 
Section 9 summarizes the obtained results and contains some suggestions regarding further 
developments. 
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2 The generalized method of characteristics 



The generalized method of characteristics has been designed for the purpose of solving quasi- 
linear hyperbolic systems of first order PDEs in many dimensions. This approach enables us 
to construct and investigate Riemann waves and their superpositions (i.e. Riemann k- waves), 
which are admitted by these systems. The main feature of the method is the introduction 
of new independent variables (called Riemann invariants) which remain constant on certain 
hyperplanes perpendicular to wave vectors associated with the initial system. This results 
in a reduction of the dimensionality of the problem. A number of attempts to generalize the 
Riemann invariants method and its various applications can be found in the recent literature 
of the subject (see e.g. [T2] - [15] , [SI], [38J and references therein). 

At this point, we summarize the version of the GMC for constructing k-wave solutions 
developed progressively in [6], [15] . [19], [39], |40j. Let us consider a quasilinear hyperbolic 
system of / first order PDEs 

A^(u)uf = 0, /i = l,...,Z, a = l,...,q, i = l,...,p, (2.1) 

in p independent variables x = (x 1 ,...,x p ) G X C W and q dependent variables u = 
(u 1 , . . . ,u q ) G U C M. q . The term uf denotes the first order partial derivative of u a with 
respect to x\ i.e. it" = du a /dx l . Here we adopt the summation convention over the repeated 
lower and upper indices, except in the cases in which one index is taken in brackets. The 
system is properly determined if I = q. All considerations have local character, that is, 
it suffices to look for solutions defined in a neighborhood of x = 0. The main steps in 
constructing k-wave solutions can be presented as follows. 

1. Find the real-valued functions X A = (X A , . . . , A^) G X and 7^ = (7^, • • • ,7^) G U by 
solving the wave relation associated with the initial system (12.11) . 

{A^{u)Xf = 0, A=l,...,k<p. (2.2) 

Thus we require that the condition 

rank {A^ a (u)Xf) < min (I, q) (2.3) 

holds. We assume here the generic case in which the rank does not vary on some open subset 
Q C U. This step is completely algebraic. 

2. Let us assume that we have found k linearly independent functions X A and 7^ which are 
C 1 in Q. We postulate a form of solution u(x) of the initial system (12.1 1) such that all first 
order derivatives of u with respect to x % are decomposable in the following way 

r)n a k 

-g-(x) = J2h A (x)r A (u)Xf(u) (2.4) 

on some open domain B C X X U. Here, h A (x) are arbitrary functions of x. This step means 
that the original system 02. ip is subjected to the first order differential constraints (12.4ft . 
Thus we have to solve an overdetermined system composed of (12.11) and (12. 4p . 
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The condition ( 12 Ah . crucial to the GMC, determines the class of solutions, called Riemann 
k-waves, resulting from superposition of k simple waves. 

3. Before proceeding further, we should verify whether the conditions on the vector functions 
X A and 7 a, which are necessary and sufficient for the existence of solutions of the system 
composed of (12. ip and (I2.4p . are satisfied. These conditions, in accordance with the Cartan 
theory of systems in involution [7], take the form 

i) [ja,1b] e span{7 A ,7 B }, ii) C JB X A E spanjA" 4 , X B }, A ^ B = 1, . . . , k, (2.5) 

where £ 7B denotes the Lie derivative along the vector field 7# and the bracket [7^, 7b] denotes 
the commutator of the vector fields ■ja,1b- 

4. Given that the conditions (12.51) are satisfied, we can choose, due to the homogeneity of 
the wave relation (12.21) . a holonomic system for the vector fields {71, . . . , 7^}, by requiring a 
proper length for each vector 7^4 such that 

[7a,7b] = 0. (2.6) 

Conditions (12.61) determine a k-dimensional submanifold S C U which can be obtained by 
solving the system of PDEs 

f)n a 

7 >V ..,««), A = l,...,k. (2.7) 



dr A 



The solution of (TJ^D 



u = f(r\...,r k ) (2.8) 

gives the explicit parametrization of the submanifold S immersed in the space U . 

5. Next we consider the functions f*(X A ), that is, the functions X A (u) pulled back to the 
submanifold S C U. The \ A {u) become functions of the parameters (r 1 , . . . ,r k ) on S. For 
simplicity of notation, we denote f*(X A ) by X A (r 1 , . . . ,r k ). 

6. Restricting the equations (12.41) and (12.51) to the submanifold S and using the linear 
independence of the vectors {71, . . . , 7^}, we obtain 

h A (x)X A (r\...,r k ), (2.9) 



dx l 
d\ A 
dr B 



a A (r\ r k )X A + (3 A (r\ r k )X B , A^B = l,...,k, (2.10) 



for some real-valued functions and (3g : S — > M. Here we do not use the summation 
convention. According to the Cartan theory of systems in involution, the conditions ( 12.71) . 
( 12. 9p and ( 12.101) ensure that the set of solutions of the initial system of PDEs ( 12.11) subjected 
to the differential constraints (12. 4p depends on k arbitrary functions of one variable. 
7. Next, we look for the most general class of solutions of the linear system of equations 
( 12.101) for X A as functions of r 1 , . . . , r k . We can perform this analysis by using, for example, 
the Monge-Darboux method 
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8. From the general solution of ( 12.101) for the functions \ A , the solution of the system ( 12. 9ft 
can be derived in the implicit form 

Af(r\...,r fe y = ^ A (r\...,r k ), A = l,...,k, (2.11) 

where ip A : R fc — > R are some functionally independent differentiable functions of k variables 
r 1 , . . . , r k such that 

() '' [ a A {r\...,r k )^ + (3 A (r\ . . . : r k )^ B \ A^B. (2.12) 



dr B 

Note that the solutions r A (x) of (12.111) are constant on (p-k)-dimensional hyperplanes per- 
pendicular to the wave vectors \ A . 

9. Finally, the k-wave solution of (12 .ip is obtained from the explicit parametrization (12. 8p 
of the submanifold S C U in terms of the parameters r 1 , . . . ,r k , which are now implicitly 
defined as functions of x , . . . , x p by the solutions of the system (12.111) in the space X . 
If the set of implicitly defined relations between the variables u a , x l and (r 1 , . . . ,r k ), 

u a = f a (r\...,r k ), Af(r\...,r fe y = ^ A (r\ . . . , r fc ), (2.13) 

can be solved in such a way that r and w a can be given as graphs over an open subset 
T> C X, then the functions u a = f a (r l (x), . . . , r k (x)) constitute a k-wave solution of the 
quasilinear hyperbolic system ( 12. ip . The scalar functions r A (x) are called the Riemann 
invariants. For p = 2 they coincide with the classical Riemann invariants as they have been 
usually introduced in the literature of the subject (see e.g. [TO] . [26] . [4"T] . [42J ) . 

Finally, let us comment on the Cauchy problem for Riemann k-waves (for a detailed 
discussion see e.g. [19].[30].[4T)]). 

Let us consider q functions u 1 (a;), . . . ,u q (x) which take some prescribed values Uq(x) = 
(ul(x), . . . ,Mq(x)) on the hyperplane Pel" defined by t = 0. Here, we use the notation 
x = (t, x) G X C lR n+1 . It was shown [H] that for < t < T the initial value problem for 
the system (12.11) has locally exactly one solution in the form of a Riemann k-wave defined 
implicitly by relations (12. 8p . (12. lip and (12.121) if the C 2 function uq{x) satisfies the following 
conditions. 

i) Uq(x) is sufficiently small that there exists a time interval [0,T] in which the gradient 
catastrophe for a solution u(x) of (12. ip does not occur. 

ii) uq(x) is decomposable according to conditions (12.41) . that is 



du a 
dxi 



x 



J2^ A (0,x)rA(Mmt(M^) (2-14) 



A=l 

on some open domain £ C P x U . 



3 Conditional Symmetries and Riemann Invariants 

Until now, the only way to approach the problem of superposition of many Riemann waves 
in multi-dimensional space was through the GMC This method, like all other techniques of 
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solving PDEs, has its limitations. They have motivated the authors to search for the means 
of constructing larger classes of multiple wave solutions expressible in terms of Riemann 
invariants. The natural way to do it is to look at these solutions from the point of view 
of group invariance properties. The feasibility and advantages of such an approach were 
demonstrated for certain fluid dynamic equations in [23], [23]. We have been particularly 
interested in the construction of nonlinear superpositions of elementary solutions (i.e. rank-1 
solutions) of (12.11) . and the preliminary analysis indicated that the method of conditional 
symmetry is an especially useful tool for this purpose. 

We use the term "conditional symmetry" here as introduced by P.J. Olver and P. Rosenau 
[33] . It evolved from the notion of "nonclassical symmetry" which had originated from the 
work of G. Bluman and J. Cole [2J and was developed by several authors (D. Levi and P. 
Winternitz [2_S] and Fushchych [TB] among others). For a review of this subject see e.g. [3J, 
[8], [35J and references therein. 

The method of conditional symmetry consists in supplementing the original system of 
PDEs with first order differential constraints for which a symmetry criterion of the given 
system of PDEs is identically satisfied. Under certain circumstances this augmented system 
of PDEs admits a larger class of Lie symmetries than the original system of PDEs. For our 
purpose we adapt here the version of CSM developed in [22], [23] . 

We now reformulate the task of constructing rank-k solutions expressible in terms of 
Riemann invariants in the language of the group theoretical approach. Let us consider the 
nondegenerate system (12 .ip in its matrix form 

A 1 {u)u l + ... + A p (u)u p = 0, (3.1) 

where A 1 , . . . , A p are I by q real- valued matrix functions of u. If we set / = q, p = n + 1 (we 
denote the independent variables by t = x ^ 1 , . . . ,x n ) and A is the identity matrix, then 
the system has the evolutionary form 

rt 

u t + J2A i (u)u j = 0. (3.2) 

i=i 

For a fixed set of k linearly independent real- valued wave vectors 

X A (u) = (Xf(u),...,X A (u)), A=l,...,k<p, 

with 

ker (XfA 1 ) ± 0, (3.3) 
we define the real-valued functions r A : X x U — > R such that 

r A {x,u) = Xf{u)x\ A = l,...,k. (3.4) 

These functions are Riemann invariants associated with the wave vectors X A , as introduced 
in the previous section. 
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We postulate the form of solution of (13. ip defined implicitly by the following set of relations 
between the variables u a , x % and r A 

u = j{r 1 {x,u), . . . ,r k (x,u)), r A (x,u) = X A (u)x\ A=l,...,k. (3.5) 

Equations (13. 5p determine a unique function u(x) on a neighborhood of x = for any 
/ : R k — > R q . The Jacobi matrix of equations (13. 5p can be presented as 

^(<H(z,-f .ly'fAeR- (3.6) 



or equivalently as 



du = ^(x k -£.-?f)~\ew*>, (3.7) 



dr \ du dr 



where 



df fdf 



dr \ dr A 
dr ( dr A \ ( d\i 



e R qxk , A = (Af) e R kxp , (3.8) 
^-x l ^eR kxq } r = (r\...,r k ) eR k , (3.9) 



du \du a J \du c 

and X q and X\~ are the q by q and k by k identity matrices respectively. Applying the implicit 
function theorem, we obtain the following conditions ensuring that r A and u a are expressible 
as graphs over some open subset T> of R p , 

*°> (3 ' 10) 

or 

det ( r '-^"l)^°- (3 - u ) 

The inverse matrix in (I3.6P (or in (13 ,7p ) is well defined, since 

dr 

= at x = 0. (3.12) 

du 

In our further considerations we assume that the conditions ( 13.101) or (13.111) are fulfilled, 
whenever applicable. 

The postulated solution ( 13. 5p is a rank-k solution, since the Jacobi matrix of u(x) has 
a rank equal to k. Its image is a k-dimensional submanifold Sk in the first jet space J 1 = 
J\X x U). 

For a fixed set of k linearly independent wave vectors {A 1 , . . . , \ k } we define another set 
of (p — k) linearly independent vectors 

Uu) = (e a (u),---,e a (u)) T , a = l,...,p-k, (3.13) 
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satisfying the orthogonality conditions 



0, A=l k. 



(3.14) 



Then, due to (13.61) (or (13. 7p ). the graph of the solution T = {(x, u(x))} is invariant under the 
family of the first order differential operators 

d 



Xa = e a (u) 



l,...,p-k, 



(3.15) 



defined on X x U. Note that the vector fields X a do not include vectors tangent to the 
direction u. So, the vectors fields X a form an Abelian distribution on X x U, i.e. 



[X a ,X b ]=0, a^b=l,...,p-k. 



(3.16) 



Conversely, if u(x) is a q-component function defined on a neighborhood of x = such 
that the graph T = {(x,u(x))} is invariant under a set of {p — k) vector fields X a with 
properties (13.141) . then u(x) is a solution of equations (13. 5p . for some /. This is so, because 



the set {r 1 , 



u q } constitutes a complete set of invariants of the Abelian algebra 



of the vector fields (I3.15p . This geometrically characterizes the solutions u(x) of the equations 

(EES). 

The group-invariant solutions of the system (13.11) consist of those functions u = f(x) 
which satisfy both the initial system (13.11) and a set of first order differential constraints 



C< = 0, a = l,...,p-k, 



a 



q, 



(3.17) 



ensuring that the characteristics of the vectors fields X a are equal to zero. 

Note that, in general, the conditions (13.171) are weaker than the DCs (12.41) required by the 
GMC, since the latter are submitted to the algebraic condition (I2.2p . Indeed, (13.171) implies 



where are real- valued matrix functions on the first jet space J 1 = J X (X x U), 



or 



df dr \ 



-i 



dr B 



' dr du ) 

dr df 



d_f_ 

dr A 



Ik 



du dr 



B 



ogxfc 



DlJXfc 



(3.18) 



(3.19) 



(3.20) 



which do not necessarily satisfy the wave relation (12. 2p . This fact results in easing up the 
restrictions on initial data at t = 0, thus we are able to consider more diverse configurations 
of waves involved in a superposition than in the GMC case. 

We now proceed to solve the overdetermined system composed of the initial system (13. ip 
and the DCs (I3TT7|) 



AZ(u)u? = 0, = 0, n 



,1, 



l,...,p-k. 



(3.21) 
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Substituting (ESJ) (or (1577I) ) into (ED) yields 



or 



9/ / <9r df 



-i 



1,| ^^"^-^J A|=0 ' /, = 1 (3 ' 23) 

where .4 1 , . . . , A 1 are p by g matrix functions of w (i.e. ^ = (A^M) e w xq , n = i,...,l). 
For the given system of equations (12. II) . the matrices A^ are known functions of u and 
equations (13.221) (or (13.231) ) constitute conditions on functions f a (r) and X A {u) (or, by virtue 
of (13.141) . on £ a (w)). It is convenient from a computational point of view to split x % into x %A 
and x %a and to choose a basis for the wave vectors X A such that 

\ A = dx iA + X A dx ia , A=l,...,Jfe, (3.24) 

where (lA^a) is a permutation of (1, . . . ,p). Hence, expression ( 13.91) becomes 

dr A d\ A . 

^- = (3.25) 

du a du a v ' 

Substituting fl3~25|) into fl3~22l) (or fl3T23|) ) yields 

tr^^-Q^-^Aj =0, /. = !,...,/, (3.26) 



or 



where 



tr (^A^(l k -K a x ia ) ^ =0, /i = l,...,/, (3.27) 



Q« = ^tfo e M 9X<3 , K a = ^-^ e M fexfc , % = I ) e IR 9 , (3.28) 



for %a fixed and i a = l,...,p — 1. Note that the functions r A and x %a are functionally 
independent in a neighborhood of x = 0. The matrix functions A 11 , df/dr, Q a and K a 
depend on r only. Hence, equations (I3.26P (or (13.271) ) have to be satisfied for any value of 
coordinates x %a . This requirement leads to some constraints on these matrix functions. 

According to the Cayley-Hamilton theorem, for any n by n invertible matrix M, the 
expression (M _1 det M) is a polynomial in M of order (n— 1). Thus, using the tracelessness of 
the expression Af 1 (T q — Q a x la ) (df /dr)X, we can replace equations (I3.26P by the following 

tr (A^Q^Xj = 0, where Q = adj(J g - Q a x ia ) G R qxq . (3.29) 
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Here adjM denotes the classical adjoint of the matrix M. Note that Q is a polynomial of 
order (q — 1) in x la . Taking (I3.29p and all its partial derivatives with respect to x ta (with r 
fixed at x = 0), we obtain the following conditions for the matrix functions f(r) and A(/) 

tr LA^x) = 0, fi = l,...,l, (3.30) 

tr^g (ai ...g as) ^AJ =0, (3.31) 

where s = 1, . . . , q — 1 and ( Cb\ , . . . , CL S ) denotes the symmetrization over all indices in the 
bracket. A similar procedure can be applied to system (13.271) to yield (I3.30P and 

tr {A^K (ai . . . K as) \^j = 0, K = adj(J fc - K a x^) G R kxk , (3.32) 

where now s = 1, . . . ,k — 1. Equations (I3.30j) represent the initial value conditions on a 
surface in the space of independent variables X, given at x %a = 0. Note that equations 
(13.311) (or (13.321) ) form the conditions required for preservation of the property (13.301) by 
flows represented by the vector fields (13.151) . Note also that, by virtue of (I3.24p . X a can be 
expressed in the form 

X a = d ia -\fd iA . (3.33) 
Substituting expressions ( 13.281) into ( 13.3 1}) or ( 13.321) and simplifying gives the unified form 



tr r ^^■■■^ A j =o ' *«={d£) e * > *= i »-' a . ( 3 - 34 ) 

where we can choose either max (s) = q — 1 or max (s) = k — 1, whichever is more convenient. 

Let us note that for k = 1 the results of the two methods, CSM and GMC, overlap. This 
is due to the fact that conditions (12.21) and (12. 7p coincide with (13.301) and conditions (12. 5p 
and (I3.34p are identically equal to zero. In this case, all rank-1 solutions correspond to single 
Riemann waves. However, for k > 2 the differences between the two approaches become 
essential and, as we demonstrate in the following examples, the CSM can provide rank-k 
solutions which are not Riemann k-waves as defined by the GMC. 

We now introduce a change of variables on ~R P x IR 9 which allows us to rectify the vector 
fields X a and simplify considerably the structure of the overdetermined system (13.211) . For 
this system, in the new coordinates, we derive the necessary and sufficient conditions for 
existence of rank-k solutions in the form (13.51) . 

Let us assume that there exists an invertible k by k subblock 

A = (A^), l<A,B<k, (3.35) 
of the matrix A G M A;><:P . Then the independent vector fields 

Xk+1 = 'dx k ~+ l ~ ( A ^^d^B ' ■ • ■ ' Xp = ~dx~P ~ ' A ^ Ap V' ( 3 ' 36 ' ) 
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have the required form ( 13. 15ft for which the orthogonality conditions ( 13.141) are satisfied. We 
introduce the functions 



x 1 = r x {x, u), . . . , x k = r k (x, u), 

™fc+l _ fc+l =p _ p -rl _ 1 T q _ q \ • J 

JU Ju . . . . • Jo Ju . Ub Uj . . . . • Uj Uj . 

as new coordinates on R p x 1R 9 space which allow us to rectify the vector fields (13.361) . So, 

we get 

d d 

Xk+1 = dx^--- Xp= dx~ P - (3 - 38) 

The p-dimensional submanifold invariant under X^+i, . . . ,X P is defined by equations of the 
form 

u = f(x\...,x% (3.39) 

for an arbitrary function / : WL k — > 1R 9 . The expression (13.391) is the general solution of the 
invariance conditions 

Ugk+i, . . . , Usp = 0. (3.40) 

The initial system (13.11) described in the new coordinates (x, u) G W x M 9 is, in general, a 
nonlinear system of first order PDEs, 

— ■ — ■ dx 1 

A i (x,u,Ux)Ui — 0, where A % = — — A> , i,j — l,...,p. (3-41) 

0X3 

That is, we have 

^ = ^—A\ ...,A h = ^—A\A k+1 = A k+1 , . . . , A p = A p . (3.42) 

ox 1 ox 1 

The Jacobi matrix in the coordinates (x, u) takes the form 

°'" ^YX G R kxp , (#) = (s) - l^l^l e M fcxfc , (3.43) 



whenever the invariance conditions ( 13.401) are satisfied. Augmenting the system ( 13.41 j) with 
the invariance conditions (I3.40p leads to a quasilinear reduced system of PDEs 

du dr \ ~ l du 



A:{ 'VA^-Si ITU) S A j= ' "= 1 --'- (3.44) 

Uj.k + 1 , . . . , UxP — 0, 

or 



A : J ^W-STSj A j= ' " =1 --'' (3.45) 
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Now we present some examples which illustrate the preceeding construction. If 

hTi ht 

is a scalar matrix, then system (13.441) is equivalent to the following quasilinear system 

B\u)ui = (3.47) 

in k independent variables x l , . . . , x k and q dependent variables u 1 , . . . , u q , where B l = XjA?. 
In the simplest case, when k — 1, the equations (I3.47P coincide with the system (13.441) . i.e. 

\ i (u)A\u)u l = ®, u 2 ,...,u p = 0, (3.48) 

with the general solution 

u(x) = fix 1 ), (3.49) 
where / : K — > M. q satisfies the first order ordinary differential equation 

K(f)AV)f = 0, (3.50) 

and we have used the following notation /' = df /dx 1 . 

If k > 2 then <p is a scalar matrix if and only if the Riemann invariants do not depend on 
the function u, 

« 3 - 5i > 

Consequently, equations f!3.25|) and (I3.5ip imply that the wave vectors A 1 , . . . , X k are constant. 
Hence, this solution represents a travelling k-wave. 

Consider now a more general situation when the matrix does not depend on variables 
x k+1 , . . . , x p , that is 

The system (13.441) is independent of x h+1 , . . . , x p if and only if 

d 2 r d 2 r 



dudx k+l ' ' dudx p 
or equivalently, due to (13.351) . if and only if 



0, (3.53) 



du du 
So, it follows that 



'' - m {k- l )™Xl, l<A<k<t<p, m,n=l,...,k. (3.54) 



— {(A-^iXT) = 0, \<A<k<i<p. (3.55) 



d_ 

du 

Thus, equations (I3.44p are independent of x k+1 , . . . , x p if there exists a A; by (p — k) constant 
matrix C such that 

(X? ) =A-C, \<A<k<i<p. (3.56) 
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In this case, (13. 44ft is a system (not necessarily a quasilinear one) in k independent variables 
x 1 , . . . , x k and q dependent variables u 1 , . . . , u q . 

Let us now proceed to define some basic notions of the conditional symmetry method in 
the context of Riemann invariants. We associate the original system (13.11) and the invariance 
conditions (13.171) with the subvarieties of the solution spaces 

Sa = {(z,u (1) ) :-4>? = 0, /i=1,...,Z} 
and 

S Q = {(x,u^):C(u)K = 0, a = l,...,q, a = 1, . . . ,p - k}, 

respectively. 

A vector field X a is called a conditional symmetry of the original system (13. ip if X a is 
tangent to § = §a H Sq, i.e. 



pr 



(i)X 



€ T (x,uW)^ ( 3 - 57 ) 



where prWX a is the first prolongation of X a defined on J l (X x U) and is given by 

pr«X a = X a -Q M/3 <<^, a = l,...,p-k, (3.58) 

and TL u (i)-j§ is the tangent space to § at some point (x,u^) E J X (X x U). 

An Abelian Lie algebra L spanned by the vector fields Xi, . . . , X p _k is called a conditional 
symmetry algebra of the original system (13 .ip if the following condition 

pi^Xai^Ui) =0, a=l,...,p — k, (3.59) 

is satisfied. 

Note that every solution of the overdetermined system (I3.2ip can be represented by its 
graph {(x,u(x))}, which is a section of J . The conditional symmetry algebra L of (13. ip 
defines locally the action of the corresponding Lie group G on J . The symmetry group G 
transforms certain solutions of (13.211) into other solutions of (13.211) . If the graph of a solution 
is preserved by G then this solution is called G-invariant. 

Assume that L, spanned by the vector fields Xi,...,X p -k, is a conditional symmetry 
algebra of the system (13.11) . A solution u = f(x) is said to be a conditionally invariant 
solution of the system (13.11) if the graph {(x,f(x))} is invariant under the vector fields 
Xi, . . . , A p _fc. 

Proposition. A nondegenerate quasilinear hyperbolic system of first order PDEs h'J.l\) in 
p independent variables and q dependent variables admits a (p — k)- dimensional conditional 
symmetry algebra L if and only if (p — k) linearly independent vector fields X\, . . . ,X p _k 
satisfy the conditions A3. 30\) and \3. 34\ ) on some neighborhood of (xq,Uo) o/S. The solution 
of A3. 1\) which are invariant under the Lie algebra L are precisely rank-k solutions of the form 
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Proof : Let us describe the vector fields X a in the new coordinates (x, u) on R p x R'. From 
( 13T38D and fl339|) it follows that 

pr^X a = X a , a = k + l,...,p. (3.60) 

Hence, the symmetry criterion for G to be the symmetry group of the overdetermined system 
(13.441) (or (13.451) ) requires that the vector fields X a of G satisfy 

X a (A) = 0, (3.61) 

whenever equations ()3.44p (or (13.451) ) hold. Thus the symmetry criterion applied to the invari- 
ance conditions (I3.40p is identically equal to zero. After applying this criterion to the system 
(13.411) in new coordinates, carrying out the differentiation and next taking into account the 
conditions (13.301) and (13.341) we obtain the equations which are identically satisfied. 

The converse is also true. The assumption that the system (13. ip be nondegenerate means 
that it is locally solvable and is of maximal rank at every point (xq,u^) G §. Therefore [36], 
the infinitesimal symmetry criterion is a necessary and sufficient condition for the existence 
of the symmetry group G of the overdetermined system (13.211) . Since the vector fields X a 
form an Abelian distribution on X x U, it follows that, as we have already shown in this 
section, conditions (I3.30p and (13.341) hold. That ends the proof, since the solutions of the 
overdetermined system (I3.2ip are invariant under the algebra L generated by (p — k) vector 
fields Xi, . . . , Xp_k- The invariants of the group G of such vector fields are provided by the 
functions {r 1 , . . . , r k , u 1 , . . . , u q }. So the general rank-k solution of (13.11) takes the form (13.51) . 

□ 

The expressions in equations (I3.30p and (13.341) lend themselves to further simplification. 
Let us recall here that any n by n matrix b is a root of the Cayley-Hamilton polynomial 

n 

det(AJ n -6) = A n -^p i (6)A n ~\ (3.62) 

i=i 

Faddeev's approach ([IT], p.87) provides a recursive method to compute the coefficients Pi(b), 
based on Newton's formulae 

n 

kp k = s k -pis fc _i - . . . -Pfe-isi, s k = tr(b k ) = ^2 k = l,...,n. (3.63) 

For example, one readily computes 

pi = tr(6), p 2 = i[tr(6 2 )-(tr(6)) 2 ], p 3 = ktr(6)) 3 - 3tr(6)tr(6 2 ) + 2tr(6 3 )], 
2 6 

p 4 = ^[6tr(6 4 ) - (tr(6)) 4 - 8tr(6)tr(6 3 ) - 3(tr(6 2 )) 2 + 6(tr(6)) 2 tr(6 2 )], (3.64) 
••• Pn = (-l) n+1 det(6). 
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According to the Cayley-Hamilton theorem one has 

n 

b n - 53^(6)6"-* = 0. (3.65) 

i=l 

Using the identity (13.651) . we can simplify the expressions (I3.30p and (13.341) to some degree, 
depending on the dimension of the matrix b. 

As an illustration we present the simplest case of a 2 by 2 matrix, which corresponds to 
rank-2 solutions for q = 2 unknown functions. In this case, the expressions (13.301) and f !3.34|) 
become 

tr^4^AJ=0, /i = l, (3.66) 
tr( A^rja^fx) =0, o = l,...,p-l. (3.67) 



dr dr 

Combining (13.661) and (I3.67P leads to the factorized form 

\ or or J \ or \ Or \ Or J J J 

Note that for any invertible 2 by 2 matrices M and N, the Cayley-Hamilton trace identity 
has the form 

MN -l 2 tr(MN) = —(N — X 2 tr(N))(M — J 2 tr(M)). (3.69) 
Using the above equation, we rewrite (13.681) in the equivalent form 



dr \dr \dr 

where the matrices M and are identified with r\ a and df /dr, respectively. Since we have 

N 2 - tr (N)N = -J 2 det (N), 
then, if det (df /dr) 0, it follows that 

/ d\ A \ 

tr(X" ( Va - l 2 tT(Va)) A) =0, Va = ^ G M 2x2 , A= 1, 2. (3.71) 

For a given system (13.11) (i.e. given functions A^), the equations ( 13. 71 ft form a bilinear system 
of l(p — 1) PDEs for 2(p — 1) functions Xf (u). Thus we have eliminated the matrix term 
df /dr in equations (13.67p . This fact greatly facilitates our task. The proposed procedure 
for constructing rank-2 conditionally invariant solutions of the system ( 13. ip consists of the 
following steps. 
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1. We first look for two linearly independent real- valued wave vectors A 1 and A 2 by solving 
the dispersion relation (13.31) associated with the initial system (13. ip . 

2. If such wave vectors do exist, we substitute them into PDEs (I3.7ip and solve this system 
for X A in terms of u. 

3. Next, we substitute the most general solutions for A 1 and A 2 into equations (13.661) and 
look for a solution u = f{r l ,r 2 ) of this system. Thus we obtain the explicit parametrization 
of the 2-dimensional submanifold S C U in terms of r 1 and r 2 . 

4. We suppose that u = /(r 1 , r 2 ) is the unique solution of PDEs f!3.66[) . Then we restrict the 
wave vectors X A to the submanifold S C U . That is, the functions \ A {u) are pulled back to 
S and become functions of the parameters r 1 and r 2 on S. We denote the function f*(X A ) 
by A^r 1 ,?- 2 ). 

5. In this parametrization we can implicitly determine the value of the Riemann invariants 
for each solution X A of (I3.7ip 

r A {x,f{r l ,r 2 )) = Xf{r\r 2 )x\ A = 1,2. (3.72) 

6. Finally, we suppose that the set of implicitly defined relations (13. 5p and (I3.72p between 
u a , r A and x % can be solved so that the functions r A and u a can be given as graphs over an 
open subset Del. Then the function 

u = f(r 1 (x),r 2 (x)) (3.73) 

is an explicit rank-2 solution of the quasilinear hyperbolic system (13.11) . The graph of this 
solution is invariant under (p — 2) linearly independent vector fields X a . 

As another illustration, let us consider the rank-3 case when q = k = 3. Then the 
condition (13.341) takes the following form 



or 



df df 



^A)=0, » = !,.. .,1, (3.74) 



where 

d\ A 



^ = I ft? I G R3X3 ' J = 1 ' 2, A=1 ' 2 ' 3 ' ( 3 - 75 ) 



and the expressions (I3.8p become 



df [df 



dr \ dr A 
We introduce the notation 



P-=Va^, Q:= V a 2 ^;. (3.76) 



Then, combining equations (13.301) and ( 13.741) . we obtain 

tv(A^[PQ-lMPQ) J rQP- X 3 tr(QP)]A ) 0. CL77I 
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If det (df/dr) ^ (otherwise the case q = 3 can be reduced to q < 2) then, similarly to 
the case q = 2, we can eliminate the term df/dr from f 1 3 . 3 j) and (I3.77p . The resulting 
expressions are still quite complicated. Nevertheless, as we show in the examples to follow, 
our procedure makes the construction of rank-3 solutions feasible. 



4 The fluid dynamics equations 

At this point, we would like to illustrate the proposed approach to constructing rank-k 
solutions with the example of the fluid dynamics equations. The fluid under consideration is 
assumed to be ideal, nonstationary, isentropic and compressible. We restrict our analysis to 
the case in which the dissipative effects, like viscosity and thermal conductivity, are negligible, 
and no external forces are considered. Under the above assumptions, the classical fluid 
dynamics model is governed by the system of equations in (3+1) dimensions of the form 

Dp + pdivu = 0, 

Du + p- 1 Vp = 0, (4.1) 
DS = 0. 

Here we have used the following notation : p, p and 5* are the density, pressure and entropy 
of the fluid, respectively, u = (w 1 , w 2 , w 3 ) is the vector field of the fluid velocity and D is the 
convective derivative 

D = j t + {u-V). (4.2) 

Equations (14. ip form a quasilinear hyperbolic homogeneous system of five equations in five un- 
known functions (p, p, u) G M 5 . The independent variables are denoted by (x M ) = {t, x, y, z) G 
X C M 4 , pL = 0,1,2,3. According to l33j,[37j this system can be reduced to a hyperbolic 
system of four equations in four unknowns u = (u^) = (a,u) G U C M 4 describing an isen- 
tropic ideal flow, when the sound velocity a is assumed to be a function of the density p 
only. In this case the state equation of the media p = f(p, S) is subjected to the differential 
constraints 

Vp=a 2 (p)Vp, d In (ap- 1/K ) = 0, (4.3) 

where a 2 (p) = df /dp, k = 2(7 — 1) _1 and 7 is the adiabatic exponent of the fluid. Under 
the assumptions (I4.3p . the fluid dynamics model (14. ip becomes 

Da + K~ l adhru = 0, ^ 
Du + /taVa = 0. 

The system of equations (14.41) can be written in the equivalent matrix evolutionary form 
(13. 2p . Here n = 3 and the 4 by 4 matrix functions A 1 , A 2 and A s take the form 



A 



/ 




8i\K 1 a 


5i2& 1 a 








SuKa 














5 i2 Ka 





v? 







V 


5 i3 Ka 








a' 


J 



1,2,3, (4.5) 
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where Sij = 1 if i = j and otherwise. The largest Lie point symmetry algebra of these 
equations has been already investigated in [20]. It constitutes a Galilean similitude algebra 
generated by 12 differential operators 



Pfi — d x v, Jk — £kij(x l d x j + u l d u j), Ki — td x i + d u i, i — 1, 2, 3, 
F = td t + x i d x i, G = -tdt + ada + rfdui. 

In the particular case when the adiabatic exponent is 7 = 5/3, this algebra is generated 
by 13 infinitesimal differential operators, namely the 12 operators (14.61) and a projective 
transformation 

C = t{td t + x% - ad a ) + (ar* - zV)^. (4.7) 

Note that the algebras generated by (14. 6 j) and by (14.61) with (14.71) are fibre preserving. The 
classification of the subalgebras of these algebras into conjugacy classes is presented in [20] . 
Large classes of solutions of the system (I4.4p . invariant and partially-invariant (with the 
defect structure 5 = 1), have been obtained in [2"T] . 

The wave vector A can be written in the form (A ,A), where A = (Ai,A 2 ,A 3 ) denotes a 
direction of wave propagation and the eigenvalue A is a phase velocity of a considered wave. 
The dispersion relation for the isentropic equations (14.41) takes the form 

det (A Q (u)Zi + Xi(u)A\u)) = (A + u ■ A) 2 [(A + u ■ A) 2 - a 2 A 2 ] =0. (4.8) 

Solving the dispersion equation (|4.8|) . we obtain two types of wave vectors, namely the 
potential and rotational wave vectors 

i) X E = (ea + u ■ e, — e), e = ±1, 

s ^ - ^ \ 2 ( 4 - 9 ) 

ii) A = ([u, e, m], — e x m), |e| =1, 

where e and m are unit and arbitrary vectors, respectively. Here, the equation (I4.9t i) has 
a multiplicity of 2. The quantity [u, e, fn] denotes the determinant of the matrix based on 
these vectors, i.e. [u, e, m] = det(w, e, m). Several classes of k-wave solutions of the isentropic 
system (I4.4p . obtained via the GMC, are known [5], [39]. Applying the CSM to this system 
allows us to compare the effectiveness of the two approaches. 

5 Rank-1 solutions of the fluid dynamics equations 

Analyzing the rank-1 solutions associated with the wave vectors X E and X s we consider 
separately two cases. 

In the first case, the potential wave vectors are the nonzero multiples of 

X E = (ea + e ■ u, -e), |ef = 1, e = ±1. 

The corresponding vector fields Xi and Riemann invariant r(x,u) are given by 

d d 

Xi = —(a + e ■ u)~ l ei— + — — , i = 1,2. 3, r(x, u) = (a + e ■ u)t — e • x. 

at ox 1 
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We can now consider rank-1 potential solutions, invariant under the vector fields {Xi, X2, X%}. 
The change of coordinates 

t = t, x = r(x,u), x = x , x = x , a = a, u — u , u = u , u = u , (5.1) 

on IR 4 x R 4 transforms the fluid dynamics equations (14.41) into the system 

^-,-i e .^! ^L- Ke .^L i- 1 2 3 (5 2) 

OX 1 OX OX OX 

with the invariance conditions 

Of = o S 3 =0, uf = u° =0, j = 2, 3, a = 1, 2, 3. 

If the unit vector e has the form e = (cosw 1 sinw 2 , cosw 1 cosm 2 , sin-u 1 ), then the general 
rank-1 solution is given by 

a(i, x) = k^x 1 + a , u 1 ^, x) = - In |Ccosw 2 |, CeM 
w 2 (t, x) = m 2 (x 1 ), w 3 (t, x) = — / tan(ln |C cos s|) cos s <is. 

In particular, if e is a constant unit vector, then we can integrate (I5.2p and the solution is 
defined implicitly by the equations 

a(i, x) = aix 1 ), u l (t, x) = ne^ix 1 ) + Cj, Cj G K, i = 1, 2, 3. 

If we choose a = A1X 1 , A\ G M and Cj = 0, then the explicit invariant solution has the form 

a(t,x) = [Ai(l + K)t - l] -1 Aie- x, u(t,x) = [Ai(l + n)t - l^nA^e ■ x)e. (5.3) 

Note that if the characteristics of one family associated with the eigenvalue Xq = a + e ■ u 
intersect, then we can choose a particular value of time interval [toj^], where T = (^4i(l + 
k)) _1 , in order to exclude the possibility of a gradient catastrophe. Hence, if the initial data 
are sufficiently small at t = t , then the solution (15.31) remains a rank-1 solution for the time 
t G [to,T), and no discontinuities (e.g. shock waves) can appear. 
In the second case, we fix a rotational wave vector 

X s = ([u, e, m], — e x m), |e| 2 = 1, 

and the corresponding vector fields (13.361) are given by 

d d 

Xi = [u,e,m}~\ex m)i— + —, i= 1,2,3. 

Hence, the Riemann invariant associated with X s has the form 

r(x, u) = [u, e, m]t — [x, e, m\. (5.4) 
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After substituting (15.41) into (15. ip . the change of coordinates transforms the initial system 
(I4.4p into the overdetermined system composed of the following equations 



= 0, (eiirij - ejrrii)— = 0, i^j = 1,2,3, (5.5) 

and the invariance conditions 

at = a& =0, uf = u% = 0, j = 2, 3, a = 1, 2, 3. (5.6) 

Hence, the sound velocity is constant (a = a^). If e and fh are constant vectors such that 
{t\m 2 — e 2 mi) 7^ 0, then we can integrate the system composed of (15.51) and (15.61) . The 
explicit solution is given by 

a(t,x) = a , u (t, x) = u (x ), u 2 (i, x) = u 2 {x l ), 

w 3 (t, x) = (e 1 m 2 - e 2 m 1 )~ l [C - (e 2 m 3 - e 3 m 2 )u 1 - (e 3 TOi - eim 3 )u 2 ] , 

where u 1 and m 2 are arbitrary functions of the Riemann invariant, which takes the form 

r(x, u) — Ct — [x, e, m]. (5.8) 

As expected, this result coincides with the solution obtained through the GMC [39]. The 
presence of arbitrary functions in the obtained solution allows us to find bounded solutions 
valid for all time t > 0. For example, the bounded bump-type solution u % = sech.(A^r l ), 
i = 1,2, contains no discontinuities. 



du ^ _ 
ox 1 



6 Rank-2 solutions 

The construction approach outlined in Section 3 has been applied to the isentropic flow 
equations (14. 4p in order to obtain rank-2 and rank-3 solutions (the latter are presented in the 
next section). In the case of rank-2 solutions, in order to facilitate computations, we assume 
that the directions of wave propagation \ A are constant, but not their phase velocities A^. 

After considering all possible combinations of the potential and rotational wave vectors 
(X Ei and X Si , respectively, i = 1,2,3) we found eight cases compatible with the conditions 
(13.301) and (13.341) . leading to eight different classes of solutions. These solutions, in their 
general form, possess some degree of freedom, that is, depend on one or two arbitrary func- 
tions of one or two variables (Riemann invariants), depending on the case. This arbitrariness 
allows us to change the geometrical properties of the governed fluid flow in such a way as to 
exclude the presence of singularities. This fact is of a special significance here since, as is well 
known [1], in most cases, even for arbitrary smooth and sufficiently small initial data 

at t — to, the magnitude of the first derivatives of Riemann invariants becomes unbounded in 
some finite time T; thus, solutions expressible in terms of Riemann invariants usually admit 
a gradient catastrophe. Nevertheless, we have been able to demonstrate that it is still pos- 
sible in these cases to construct bounded solutions and, in particular, soliton-like solutions, 
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through the proper selection of the arbitrary function(s) appearing in the general solution. 
To this purpose we submit this arbitrary function(s), say v, to the differential constraint in 
the form of the nonlinear Klein-Gordon equation 

v tt -Av = cv 5 , ceR (6.1) 

which is known to possess rich families of soliton-type solutions (see e.g. [I], [13]). Equation 
(16. ip can be reduced to a second order ODE for v as a function of a Riemann invariant and 
can very often be explicitly integrated. The analysis of the singularity structure of these 
ODEs allows us to select soliton-like solutions for v which, in turn, in many cases, lead to 
the same type of rank-2 and rank-3 solutions of the system (14. 4p . Among them we have 
various types of algebraic soliton-like solutions (admitting no singularity other than poles), 
kinks, bumps and doubly periodic solutions which are expressed in terms of Jacobi elliptic 
functions. 

Below we list the obtained rank-2 solutions. Some of the general solutions found (both 
rank-2 and rank-3 in the next section) coincide with the ones obtained earlier by means of 
the GMC. Nevertheless, we list them all since we derive from them the particular bounded 
solutions, which, to our knowledge, are all new. 

For convenience, we denote by {E^E^ EiSj, SiSj, EiEjE k , . . . ,i,j,k — 1, 2, 3) the solutions 
which result from nonlinear superpositions of rank-1 solutions associated with given wave 
vectors \ Ei or \ Si . The sign (+ or — ) coincides with the value of e = ±1 in equation (14.91) . 

Case (EiE 2 ). We first discuss the superposition of two potential rank-1 solutions Ei for 
which the wave vectors have the form 

\ E t = (ea + e i -u,-^), \e l \ 2 = 1, i = 1, 2, e = ±1. (6.2) 

We assume that the wave vectors A El± and X E2± are linearly independent. The corresponding 
vector fields (13.361) are given by 

X —-^l—^^l— x-— / fi ox 

1 dx 2 f3 1 dt p x dx 1 ' 2 dx 3 p x dt fix dx 1 ' 1 ' ' 



with 



Pi = e 2 (a + e 1 • u) — e\(a + e 2 • u), i = 1, 2, 3, 



a J = e \ e ) ~ e i e i> 3 = 2 > 3 > e = !• 
The nonscattering rank-2 potential solution (EfE^) has the form 

a = ai(r r ) + a 2 (r 2 ), u = K(ai(r 1 )e 1 + a2(r 2 )e 2 ), (6.5) 
where a\ and a 2 are arbitrary functions of the Riemann invariants 



(6.4) 



r 1 ^,^) = (1 + K)ai{r l )t - e 1 • x, \e l \ 2 = 1 
r 2 (x,u) = (1 + K)a 2 (r 2 )t - e 2 ■ x, \e 2 \ 2 = 1 



(6.6) 
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respectively, and the wave vectors e 1 and e 2 have to satisfy the condition 

e 1 • e 2 + /sT 1 = 0. (6.7) 
Equation (16.71) holds if and only if the angle <p between these vectors is 

cosy>= -(I-7). (6.8) 

This solution represents a Riemann double wave. Here, the rank-1 solutions Ef , i = 1,2, 
do not influence each other (they superpose linearly). This result coincides with the one 
obtained earlier by means of the GMC [59] . 

i) In the particular case when cij(r l ) = —A^r 1 , A± £ R, i — 1,2, the solution (16.51) takes 
the explicit form 

Aie 1 ■ x A 2 e 2 ■ x 

+ 



1 + (1 + K )Ait 1 + (1 + K)A 2 t 

KA\e • x _n KA 2 e ■ x _^ 

u = e -\- c 

l + + \ + {l + K)A 2 t 

which admits the gradient catastrophe at the time t = min (A7 1 (l + K )- 1 ),i = l,2. Hence, 
some discontinuities can occur e.g. shock waves which correspond to the formation of a 
condensation jump from the compression waves related to X Ei . 

ii) The following bounded solution can be obtained using the DC (16.11) 



a = Ar l (l + B^) 2 )- 1 / 2 , At, B t £ R, B, t > 0, 

i=l 
"2 

^^(1 +B i (r < ) 2 )" 1/2 e i 



U = K 



(6.10) 



j=i 

where the Riemann invariants are given by 

r i =[(l + K)A i r i (l + B i (r i ) 2 )- 1/2 ]t-e i -x, % = 1,2. (6.11) 

The result (16.101) represents an algebraic kink-type solution which is bounded for t > while 
each r l possesses a discontinuity at time T = (Aj(l + k)) . 

Case (E1S2). In the mixed case (EfS 2 ), we consider the superposition of the rank-1 
potential solution with the rank-1 rotational solution S 2 associated respectively with the 
wave vectors 



\ E t = (a + e 1 ■ m, -e 1 ), 

X s2 = ([u,e 2 ,m 2 ],-e 2 x m 2 ), |e*| 2 = 1, z = l,2. 
The vector fields (I3.36P corresponding to the wave vectors (16.121) are 



(6.12) 



d a 2 d (3 2 d d cr 3 d (3 3 d 

1 = Ih 2 ~ ft dt ~ J x dx 1 ' 2= dx^~ 'Jidi~ '^hdx 1, ^' ' 
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where 

Pi = — (e 2 x rh 2 )i(a + e 1 • u) + e^u, e 2 , rh 2 }, i — 1, 2, 3, 



e}(e 2 x m 2 )^ + ej(e 2 x m 2 ) 1; j = 2, 3. 



(6.14) 



The invariant nonscattering rank-2 solution (.E^S^) has the form 

a = cti(r 1 ) + a , m = fta^r 1 )*? 1 + M 2 (r 2 ), (6.15) 

where 

[w 2 , e 2 , m 2 ] = C 2 , |ef = 1, 2 = 1, 2, (6.16) 

and ai and u\ are any different iable functions of r 1 and r 2 , respectively, and the relation 
u 2 (r 2 ) = Ciul(r 2 ) holds. Here, ao,Ci,C 2 G R and m 2 is an arbitrary constant vector. The 
wave vector X s2 takes the form 

A* 2 = (C 2 , -(ejej + C x (l - (e}) 2 )), -e&J - C^ 1 ), (1 - (e^) 2 + C x e\e\)). (6.17) 

From fl6U5D . fl6TT6|) and (16371) . we get 

[e\e 2 ,m 2 ] = 0, (6.18) 

so the vector \ E i = — e 1 is orthogonal to A 52 = — e 2 x m 2 . Hence, the Riemann invariants 
are given by 

r 1 = ((1 + n) ai (r l ) + C 2 (C l e\ - e\)- l )t - e 1 • x, (6.19) 
r 2 = C 2 t - {e\e\ + C 1 {\ - (ejf)^ 1 - e *(eJ - C ie \)x 2 + (1 - (4) 2 + deje^x 3 . 

This solution represents a Riemann double wave. 

i) An explicit form of the solution (16.151) can be found when e 1 = e 2 = (cosy?, sin 0) 
and m 2 = (siacp, — cosy?, C\ sin<£>), Ci G R, and we choose ai(r 1 ) = Air 1 , A\ G R. The 
Riemann invariants are now given by 

j C 2 t + ar'-C]. cos 2 <£> + x 2 C*i cos^sin^ 
= (Cicos^^^l + K)t- 1) ' (6.20) 
r 2 " ' 



2 - C 2 t — Ci^ 1 sin 2 ip + x 2 C\ sin cos ip + x 3 , 



and the solution becomes 



Cot + C\X X cos 2 + CiX 2 cos ip simp s „ X/ 2 \ 
(Ci cos<^J(Ai(l + k)£ — 1) 
i K,Ai(C 2 t + Cix 1 cos 2 ip + C\X 2 siiaip cos ip) x . 2 

" = c^d + ^-i) + "> {r >' 



u 



2 



kAi 



(6.21) 



(C 2 tan ipt + Cix 1 sin cos ip + Cix 2 sin 2 ip) 



Ci 

C 2 i . o\ 

u 2 (r ) cot v?, 



C\ sin y> cos <^ 
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where u 2 (r 2 ) is an arbitrary function of r 2 . Note that a and u 1 admit the gradient catastrophe 
at the time T = (Ax(l + k)) _1 . 

ii) Another interesting case of a conditionally invariant solution occurs when we impose 
condition (16. ip on the functions a\ and u\. Then the solution is bounded and represents a 
solitary double wave of the type (E±S 2 ) 

a = Ai{l + B 1 {r 1 - l) 2 )" 1 / 2 + a , A 1 ,B 1 ,C 1 eR, B 1 >0, 

u = K A l (l + B l (r 1 - l) 2 )- l ' 2 e l + (^(r 2 ), E 2 u\(r 2 ) + F 2 , C lU \(r 2 )) T , !<> ' 

where 

u 1 2 {r 2 ) = A 2 (l + B 2 coshD 2 (r 2 -l))- 1 / 2 , A 2 ,B 2 ,D 2 eR, B 2 >0, (6.23) 
E 2 = -(eKC.el - el))- 1 ^ + e})(C ie } - ej), F 2 = C 2 (4(C ie } - e^))" 1 . 

The Riemann invariants take the form 

r 1 = (1 + «)(A X (1 + ^(r 1 - I) 2 )" 1 / 2 + C 2 (C ie i - e^ 1 )* - e 1 • x, (6.24) 
r 2 = C 2 t - (ejej + d(l - (el) 2 ))^ - e*(ej - de^x 2 + (1 - (e*) 2 + Cie\e\)x* . 

The solution remains bounded even though the function r 1 admits the gradient catastrophe 
at the time T = (1 + Bxf' 2 [(1 + k)AiSi] -1 . 
Case (S1S2) : i) Let us assume that 

e 1 ^ (0,0,1), m 1 = (0,1,0), e 2 = (1,0,0), m 2 = (0,0,1). 

Then the wave vectors (14.91 i) are given by X Sl = (—u 1 , 1,0,0) and A 52 = (— w 2 , 0,1,0) and 
are linearly independent. So we are looking for rank-2 solution (SiS 2 ) invariant under the 
vector fields 

Xi = ^ + m1 7TT+ m2 ^' X * = 7T1- (6.25) 
dt dx 1 ox ox 6 

The corresponding Riemann invariants are 

r (x,u) = x — u t, r (x,u) = x — u t. (6.26) 
The change of coordinates 



t — X — X IX X — X UL X — X , 

— —1 1 — 2 2—3 3 

a = a, u — u , u — u , u = u , 



(6.27) 



transforms the system (I4.4p in this case into the equations 

cte 1 du 2 du 1 du 2 du 1 du 2 

dx 1 dx 2 ' dx 1 dx 2 dx 2 dx 1 

da da da du 1 

0, — = 0, i = 1,2,3. 



(6.28) 



dx 1 dx 2 dx 3 dx 3 
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The solution of system (16.281) has the form 

a = do, u l (t, x) = — tt-^-, u 2 (i, x) = -Jrr-i u 3 (i, x) = vPix 1 , x 2 ), (6.29) 

ox 1 ox 1 

where the function ip satisfies the homogeneous Monge-Ampere equation 

tI>&&tI>&& ~ = °> ( 6 - 30 ) 

and u 3 is an arbitrary function of two variables. Note that this solution has rank 2 but it is 
not a Riemann double wave. 

i) The proper selection of the function ip transforms the solution (I6.29p into 



a(t,x)=a , u 1 = (1 - n) f - - - — j , neZ\{l}, 

, H 3 (t,i) = M 3 (a; 1 -H 1 f,x 2 -M 2 f). 



2 2j_\ 1— n 
2T — U t N 



(6.31) 



u 2 = —n , 

x L — u t 

For n = 2, we obtain an explicit solution of the form 

a = oo, u 1 = -2-H- 2 [xH + (x 2 ) 2 ± x 2 ((x 2 ) 2 + Atx 1 ) 1 ' 2 ], 
u 2 = t~ l [x 2 ± {{x 2 ) 2 + to 1 ) 1 / 2 ], u 3 = u 3 {x l - uH, x 2 - u 2 t). 



(6.32) 



with a singularity at t = 0. 

ii) Another example worth considering is the case when fluid velocity can be decomposed 
as follows u = U\(r l ) + u 2 (r 2 ). Then we get the scattering nonsingular rank-2 solution 



■ Wj-A^r 1 )) /Ag A 2 / 

A 2 A 1 _ A 1 A 2 ( 6 - 33 ) 

u 2 = u 2 1 (r 1 )+u 2 2 (r 2 ), u 3 = u 3 (r 1 ) + 7 1 u 2 2 (r 2 ), a = o , 77 = ttt| rfrf, 

A 1 A 3 ~~ A 3 A 1 

where we introduced the notation A Si = (A ,A*), i = 1,2. The above solution is invariant 
under the vector fields 

9 a 2 d (3 2 d d a 3 d f3 3 d 

1 = dx^ ~ J x di ~ ~pl dx 1, 2= "dx^~ 'J x di~ Jidx 1, ( ' ' 

with 

at = \{\ 2 - A 1 A 2 , fa = A 2 [w, e 1 , m 1 ] - \][u, e 2 , m 2 }, i = 2, 3, j = 1, 2, 3. (6.35) 

Here u\ and wf are arbitrary functions of r 1 , w 2 , is an arbitrary function of r 2 and C» = 
[{?i(r ? ), e l , m l ], z = 1,2. The Riemann invariants take the form 

r 1 = (d + C 2 A}/A 2 )t - A 1 • f, 

r 2 = (<7 a + |d + (A 2 - u 2 ^) + (A 2 - ^(r 1 )) t - A 2 ■ x. ^ 
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Note that the Riemann invariant r 2 depends functionally on r . This means that the inter- 
acting waves influence each other and superpose nonlinearly. The result is a Riemann double 
wave. 

iii) By submitting the arbitrary functions uf, u\ and u 2 to the DC (16.11) we can construct 
the rank-2 algebraic kink-type solution of the form 

u 1 = (Xl)- 1 ^ - \\A 2 r\l + B 2 {r l f)-^ - \\A,r\l + B^f)-^} 
- (A 2 )- 1 !"^ + (A§77 + Xl)A x r 2 {l + B 1 (r 2 ) 2 )- 1 / 2 ], A, B, e R, 

i 



u 2 = A ir 2 {\ + B^r 2 ) 2 )- 1 ' 2 + A 2 r\l + B 2 {r 1 ) 2 )- 1 / 2 , B, t > 0, 



1,2,3, 



(6.37) 



u 3 = A 3 r\l + B^r 1 ) 2 )- 1 ' 2 + rjA ir 2 (l + B^r 
where the Riemann invariants are given by 

A! 

„2 



2\2\-l/2 



a — do, 



r 1 = (d + C 2 ^)t - A 1 • x, 



C 2 + d| + (a 2 - ^) A 2 r\l + B^r 1 ) 2 )^ 2 
^^ 3 r 1 (l + S 3 (r 1 ) 2 )- 1 / 2 



(6.38) 



+ A - 



AfAj 



t - X ■ x. 



Case (E1E2S3) . The nonscattering rank-2 solution (E^E^S^) invariant under the 
vector field 

d e ijk e}e 2 (e 3 x m 3 ) k d : p 23 d , /3 31 9 



A' 



<9x 3 



+ 



+ 
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dt /3i2 dx l (3i2 dx 2 



with 



e*e 2 )[u, e 3 , m 3 ] + (e 2 (e 3 x m 3 )i - e 2 (e 3 x m 3 ) j )(a + e 1 • w), 



+ (e-(e 3 x m 3 )^ - e^e 3 x m 3 )i)(a + e 2 • u), i,j = 1,2,3, 



(6.39) 



(6.40) 



has the form 



A i(( e l + ej)x 1 + (e\ + e 2 )x 2 ) 



U = Un 



u 



u 



1-Ai(1+k)< ' ' °' 

-kA x (((e\) 2 + (e 2 )V + ( e i e i + e 2 e 2 )x 2 ) - « 3 (r 3 ) 



1 - Ax(l + K)t 

e 2 (/?w 3 ( r3 )^ ~ e!^ 1 _ e^ 2 ) 



(6.41) 



1 - Ai(l + «)t 



e 2 . (— (3u\{r 3 )t — e\x l — e 2 x 2> 



1 - Ai(l + K)t 



+ 



P 2 _ 1 

c 2 21/ 3> 



where le 1 ! 2 



?2|2 



1, e 1 • e' 



-K 



-1, e J = e 2 = 0, /3 = (1 + KT 1 )/^ - e 2 ) and A 1; w 3 e K. 



The Riemann invariants are 



I3u\(r )t — e\x — e\x 
1 - A x (l + K)t 



—(3u\(r 3 )t — e\x l — e\x 2 
1 - A^l + K,)t 
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. ( r 3 = x 3 - u 3 Q t, 



(6.42) 



where u\ is an arbitrary function of r 3 . 

This solution represents a Riemann double wave. It does not admit removable singularities 
for any choice of ul(r 3 ), but the functions a, u 1 and u 2 are subject to the gradient catastrophe 
at the time T = {A x (l + k)) _1 . 

Case (E1S2S3) . The nonscattering rank-2 solution (E1S2S3) invariant under the 
vector field 

x _ d | e ljk e](e 2 x m 2 )j(e 3 x m 3 ) k d | /3 23 d | /? 3 i g , g ^ 



with 



f3ij = [(e 2 x m 2 )i(e 3 x m 3 )^ - (e 2 x m 2 )j(e 3 x m 3 )j](a + e 1 ■ u) 

+ [e](e 3 x m 3 )^ - e^e 3 x m 3 )j} [u, e 2 , m 2 ] (6.44) 
+ [e-(e 2 x m 2 ).,- -e](e 2 x m 2 )i][u,e 3 ,rh 3 }, 



is given by 



_ (CyA^ + Cg/A^t-x 1 x _ (C 2 /A 2 + C3/A 3 ) (1 - Aig - /cAxs 1 

(2 — ./ii ■ — ; : . 11 — ■ — : : . 

1-A l (l + K)t l-A x (l + K)t 

2 _ n(U\1 \3\/\2„2 , \2 3\ „,3 _ CA^&A 2 — A 3 ) . 2 9 , \2J\ 



(6.45) 



C(6A 2 - A?)(A 2 x 2 + A 2 x 3 ), u = - 2V x2 3 17 (A 2 x 2 + A 2 x 3 ). 

A 3 A X 

The Riemann invariants have the explicit form 



l-Ai(l + «)t A? 

Here, we introduced the notation A 5 ' = (Aq,A 1 ), Cj = [«i,e ! ,m ! ], i = 2, 3 and A 3 = 
(A 3 , 6A|, 6A3), ii G R . Note that a and w 1 both admit the gradient catastrophe at the 
time T = (Ai(l + k)) _1 while u 2 and u 3 are stationary. In this case the solution again has a 
form of Riemann double wave. 

Case (S1S2S3) . The rank-2 solution is invariant under the vector field 

dt ~^ dx 1 ~^ U dx 2 ~^ U dx 3 ^ 

In this case subjecting the initial system (14. 4p to the DCs (13.171) leads to the overdetermined 
system 

a = a , u + (u-V)u = 0, V5 = 0, a G R. (6.48) 
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The solution of (16.481) is divergence free if and only if 

i?=/(r 1 ,r 2 ,r 3 ), / : K 3 -> M 3 , r i = x i -u i t, i=l,2,3. (6.49) 

The Jacobi matrix Df(r) = (df a /dr l ) has to be nilpotent. In fact, the reduced system (I6.48P 
mandates that the characteristic polynomial is equal to 

det(AJ 3 + Df(r)) 

1 (6 50) 

= A 3 - A 2 tr(/° ) + -[(tr(/° )) 2 - tr(/") 2 ]A + det (/£) = A 3 . 

In order to satisfy this condition we can select the arbitrary functions f a in the following 
way f 1 = b(r 2 , r 3 ) and f 2 = f 3 = g(r 2 — r 3 ). Then we have 



Df(r) = I g s -g s I, s = r 2 - r 3 . (6.51) 
V & / 

If b r 2 ^ b r s, then r&nkDf(r) = 2, otherwise f 1 is an arbitrary function of one variable, i.e. 
f 1 = h[r 2 — r 3 ), and rank .D/(r) = 1. In the rank-2 case the solution has the form 




u\x,t) = b(x 2 - tg(x 2 - x 3 ), x 3 - tg(x 2 - x 3 )), , ( . _^ 

u 2 (x, t) = u 3 (x, t) = g{x 2 — x 3 ), a = a , ao £ m 



where b is an arbitrary function of two variables (x 2 — u 2 t) and (x 3 — u 3 t), and g is an arbitrary 
function of (x 2 — x 3 ). 

Depending on the choice of the arbitrary functions, the relations (16.521) can lead to ele- 
mentary solutions (constant, algebraic, with one or two poles, trigonometric, hyperbolic) or 
doubly periodic solutions which can be expressed in terms of the Jacobi's elliptic functions 
sn, cn and dn. To ensure that the elliptic solutions possess one real and one purely imaginary 
period and that, for real argument r l , they are contained in the interval between —1 and +1, 
the moduli k of the elliptic functions have to satisfy the condition < k 2 < 1. An example 
of such elliptic solution has been obtained by submitting the arbitrary functions b and g to 
the DC (16.11) . It has the explicit form 

u 1 = A l [l + B lS n 2 ((3(x 2 + nx 3 ) - (n + 1)[A 2 [1 + B 2 sn 2 ((3(x 2 - x 3 ), k)]~ 1/2 }, k)}' 1/2 , 
u 2 = u 3 = A 2 [l + B 2 sn 2 {(3{x 2 - x 3 ), k)}- 1 / 2 sn{f3{x 2 - x 3 ), k), (6.53) 
a = a , < k 2 < 1, Ai,Bi,PeR, B { > 0, z = l,2. 

This is a bounded solution representing a snoidal double wave. 

7 Rank-3 solutions 

Let us now present the rank-3 solutions obtained by way of the procedure analogical to the 
one described in Section 3 for the rank-2 solutions. 
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Case (E1E2E3) . The rank-3 potential solution (E^E^E^) invariant under the vector 



field 



X 



d [e\e 2 ,e 3 ] d Pi_d_ d 



dx 3 (3 3 dt ' fodx 1 ' fodx 2 ' ^ 

with /3i = (e 2 x e 3 )i(a + e 1 ■ u) + (e 1 x e 3 )i(a + e 2 ■ u) + (e 1 x e 2 )i(a + e 3 • u), takes the form 

a = a^r 1 ) + a 2 (r 2 ) + a 3 (r 3 ), u = n^a^r 1 ) + e 2 a 2 (r 2 ) + e 3 a 3 (r 3 )), (7.2) 

where the Riemann invariants are 



1 + K)a,i(r l )t - e l ■ x, e 1 ■ e J 



-K 



-1 



*i\2 



1, z^ J = l,2,3, 



(7.3) 



and a,j are arbitrary functions of r l . Note that, just as in the case E1E2, the angle (p^ 
between each pair of wave vectors e\e J , i ^ j — 1,2,3, has to satisfy the condition (16.81) . 
This nonscattering rank-3 solution coincides with the one obtained previously by the GMC 
[39]. 

i) After submitting the arbitrary functions aj to the DC ( 16. ip we obtain several bounded 
solutions. We list here two examples. 

An interesting case is the algebraic kink solution 



a = J2^ l (l + B l (r i fr 1/2 , 



U = K 



i=l 



Y J A l r i {l + B l {r l f)- 1 l 2 e l 



where the Riemann invariants are given by 

r l = [(1 + K)A t r\l + B^r 1 ) 2 )- 1 ' 2 ] t - e 1 ■ x, A i: B t G 



1,2,3. 



(7.4) 



(7.5) 



This solution evolves as a triple wave and is bounded even when the Riemann invariants r l 
admit the gradient catastrophe at the time Tj = (1 + . 

ii) Another interesting solution describes an algebraic solitary triple wave of a kink type 



-1/2 



u 



Kj2Ml + e B ' r Y 1/2 e l , Ai, Bi e 



(7.6) 



i=l 



i=l 



where the Riemann invariants are given by 

'l + K ) Ai {l + e B ^)- 1 ' 2 



t - e 1 ■ x, 



1,2,3. 



The Riemann invariants admit the gradient catastrophe at the time 

T . = _ 2 5/2 ((l + K)AiBi) -\ 



(7.7) 



{71 



but the solution remains bounded. In both cases the angle </?jj between the wave vectors e ! 
and e J is given by (16. 8p . 
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Case (E1S2S3) : In this case we have to distinguish two situations, depending on the 
choice of wave vectors X El , X s2 and A 53 . 

First we look for the rank-3 solution (E1S2S3) invariant under the vector field 

where we have assumed that the linearly independent wave vectors associated with the waves 
Ei, S2 and S3 are given by 

A s i + = (a + m 3 , 0,0,-1), A 52 = [e\u x - e*w 2 , -e 2 , 4, 0), A 5s = (-e 3 w 3 , 0, 0, e 3 ). (7.10) 

The corresponding Riemann invariants satisfy the following relations 

r 1 = ((1 + n )f(r 1 ) + a + ul)t — x 3 , r 2 = t — x 1 sing(r 2 , r 3 ) + x 2 cosg(r 2 , r 3 ), 

where r 3 obeys the evolutionary partial differential equation 

7* +W) = °' (7.11) 
The solution then takes the form 



a = k 1 f(r 1 ) + oo, u 1 = sin#(r 2 ,r 3 ,, 
M 2 = -cos5f(r 2 ,r 3 ), u 3 = f(r l ) + u 3 , a ,«o GF: 



where / is an arbitrary function of r 1 and g is an arbitrary function of r 2 and r 3 . This 
scattering rank-3 solution has been obtained earlier through the GMC [39J. 

i) If f(r l ) = A^ 1 + Bi, then the solution of (17.111) can be integrated in a closed form 



a = K 1 (A 1 r 1 + Bi) + a , u l = sing(r 2 , r 3 " 1 



2 _ „/„2 „,3\ 



A^ 1 + 5i + w 3 , A 1 ,B 1 e 



(7.13) 



and the Riemann invariants are given by 

1 ((1 + K-^B^ao + uDt-x 3 



r 



1 - (1 + K-^Arf 



r 2 = t — x 1 sing(r 2 , r 3 ) + x 2 cosg(r 2 , r 3 ), (7-14) 
r 3 = * — (^(koo - u 3 )t + x 3 - m - 5i)((l + K)Ait - k)~ k/k+1 

where \I/ is an arbitrary function of its argument and g is an arbitrary function of two variables 
r 2 and r 3 . This solution corresponds to a scattering Riemann triple wave. 

After subjecting the arbitrary functions / and g, appearing in the solution (I7.12p . to the 
DC (16.11) we get several bounded solutions. Below, we present two of them. 
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ii) A physically interesting subcase of (E1S2S3) is the solution 
a = k~ 1 A 1 [1 + + cosher 1 ))]- 1 / 2 + a , 



u = sin 



u = — cos 



A 2 (R)- 1/2 taay 
_{B 2 + tan 2 y) 1 / 2 
A 2 {R)~ 1 ' 2 tan y 



Ai, B h C x E R, Bi > 0, 



1,2, 



(7.15) 



(B 2 + tan 2 y) 1 / 2 _ 
u 3 = A 1 [1 + B x (l + cosh(CV))]- 1 / 2 + ug, 



Here we introduced the following notation R = (r 2 ) 2 + (r 3 ) 2 and y = |ln \D X R\, with £>i G 
The Riemann invariants r 1 and r 2 are 



;i + n- 1 )A 1 [l + B X {1 + cosh(C 1 r 1 ))]- 1/2 t - a; 3 , 



r 2 = t — x 1 sin 



A 2 (R)' 1/2 tany 
(B 2 + tern 2 y y/ 2 



+ X cos 



A 2 (R)~ 1/2 tany 
_(5 2 + tan 2 y y/ 2 



(7.16) 



and r 3 satisfies the linear partial differential equation 



— + (A 1 [l + B x (l + cosh(C 1 r 1 ))]-^ + u 3)_ = o. 



(7.17) 



This solution is finite everywhere except at R = 0, but has discontinuities for In = 
(2n + l)7r, fiGZ. It remains bounded even when the Riemann invariants r 1 ,?" 2 and r 3 tend 
to infinity. Physically, this solution represents nonstationary concentric waves damped by 
the factor R- 1 ' 2 . 

iii) Another solution worth mentioning has the form of an algebraic solitary wave 



a = K~ x A x \i + £1(1 + cosh Cir 1 )]- 1 / 2 + a , B X ,C X > 0, 
u l = sm(D 1 [l + e h{r2 ' r \ l/2 ), A x , B X: C x , D x e M, 
n 2 = cos( J D 1 [l + e^ 2 ^]- 1 / 2 ), 
u 3 = Ai[l + 5i(l + coshC 1 r 1 )]- 1/2 + m 3 , 

where /i is an arbitrary function of r 2 and r 3 . The Riemann invariants are 

r 1 = (1 + + B x (l + cosh(Cir 1 ))]- 1/2 t - x 3 , 

r 2 =t-x 1 smD x [l + e h ^^'\ 1 ' 2 + x 2 cosD x [l + ^(-V 3 )]-^ 



(7.18) 



(7.19) 



and r 3 satisfies the partial differential equation (17.171) . 

We now consider the case (E x S2S3) with a different selection of the wave vectors than 
assumed in f)7.10p . namely we choose 



\ E t = (a + eV + e^.-eij-e^O), |ei| 2 = 1, 
A 52 = (u 2 , 0, -1, 0), \ Ss = {-u\ 1, 0, 0). 



(7.20) 
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This leads to scattering rank-3 solution of the form 



a = K 1 f( y r 1 ) + a , u 1 = sin/(r 1 ), u 2 = — cos /(r 1 ), 
u 3 = g(r 2 cos /(r 1 ) + r 3 sin /(r 1 )), a G R, 



(7.21) 



in which g is an arbitrary function of one variable r 2 cos /(r 1 ) + r 3 sin f(r r ) . The Riemann 
invariants are 



r 1 = (k 1 f{r l ) + a )t — x 1 cos f(r r ) — x 2 sin /(r 1 ) 
r 2 = — tcos/(r 1 ) — x 2 , r 3 = — tsin/(r 1 ) + x . 



(7.22) 



This triple wave solution coincides with the one obtained through the GMC 

iv) As previously, we constructed particular solutions from (I7.2ip by requiring that the 
arbitrary function / satisfies the DC (16. ip . One of the interesting examples is a periodic 
solution 



a = k~ 1 A 1 (1 - B x cofiCir l Y l/2 + a , 

u 1 = sinAi(l - B x cosC x r l )- l/2 , A X ,B X ,C X E R, 

m 2 = -cosAi(l - fiiCosCxr 1 )- 1 / 2 , |Bi| < 1, 

u 3 = g(r 2 cos A x (l -B x cosCir 1 )" 172 + r 3 sinA^l - B x cos dr 1 )' 112 ). 
with the Riemann invariants 

r 1 = (k _1 Ai(1 - B x cosdr 1 )- 1 ' 2 + a )t - x 1 cosA x (l - B x co$C x r l )- 1/2 



(7.23) 



-x 2 sinAi(l -B x cos^r 1 )- 1 / 2 , 
~ 2 = -tcos^i(l - B x cos^r 1 )- 1 / 2 -x 2 , 
- 3 = -*sinAi(l - B x cosC x r l )- l/2 + x 1 . 



(7.24) 



This solution remains bounded even when the Riemann invariants admit a gradient catas- 
trophe. 

8 Rank-k solutions of fluid dynamics equations 

Let us now consider the isentropic flow of an ideal and compressible fluid in the case when 
the sound velocity depends on time t only. The system (14.41) in (k + 1) dimensions becomes 

u t +(u- V)« = 0, 

a t + K _1 adivu = 0, a x j = 0, j = l,...,fc, a > 0, k = 2(7 — 1)~ . 

We show that in this case our approach enables us to construct arbitrary rank solutions. 
The change of coordinates on R fc+1 x R fc+1 

i=t, x 1 = x 1 - uH, . . . ,x fc = x k - u k t, a = a, u = u<ER k , (8.2) 



31 



transforms ( 18. ID into the system 

^ = 0, ^ + K - 1 atr((X fc + tDM(x))- 1 J DM(x)) =0, — = 0, (8.3) 

where Du(x) = du/dx G R fexfc is the Jacobian matrix and x = (x l ,...,x k ) G R fc . The 
general solution of the conditions du/di = and da/dx = is 

tZ(t,x) =/(£), a(t, x) = a(t) > (8.4) 



for any functions / : R fc — > R fc and a : R — > R, respectively. Making use of (18.41) and of the 
trace identity 

9 =(lndetB) ^trfs" 1 ^ ), (8.5) 



<9r 7 V 0*, 

where B = (l k + Wf(x)) and D/(x) = f-(Z fc + tDf(x)), we obtain from (Q 

|L [In (\a(t)\ K det (l k + = 0. (8.6) 
Differentiating (18.61) with respect to x gives the condition on the flow velocity f(x) 

_ [In (det (X k + iDf(x)))} = 0- (8.7) 

Consequently, we have 

det (l k + tDf(x)) = a(x)P(i), (8.8) 

where a and f3 are arbitrary functions of their argument. Evaluating (18.81) at i — implies 
a(x) = /?(0) _1 . Therefore, 

det(l k + Wf(x)) = ^ 

and we obtain 

-^det(l k + iDf(x)) = 0. (8.9) 

Equation (18.91) holds if and only if the coefficients p n , n = 0, . . . , k — 1, of the characteristic 
polynomial of the matrix Df(x) are constant. Thus the general solution of (18. ip is 

u(t,x) = f(x), a(t) = A 1 (l+p k ^ 1 i+...+ Po i k )~ 1/K , A\ G R + , (8.10) 
with the Cauchy data 

t = 0, u{0,x) = f(x), a(0) = A 1 . (8.11) 
In the original coordinates (x, u) G R p x R 9 this rank-k solution takes the form 

u = f(x l -uH,...,x k -u k t), a(t) = A 1 (l+ Pk _ 1 t + ...+p t k )' 1/K . (8.12) 
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Note that the sound velocity a is constant if and only if the Jacobian matrix Df(x) is 
nilpotent, i.e. 

det(-XI k + Df(x)) = (-\) k . 

As an example let us consider the particular solution of (18. ip for k = 2. It is invariant 
under the vector fields 

The requirement that the coefficients p n of the characteristic polynomial (I3.64p of the Jacobi 
matrix Df(x) are constant means that 

det (D(f(x))) = B h tr(D/(z)) = 2C 1: B^de R, 

where we denote B\ = po and 2C% = p±. Solving the above conditions gives us the general 
rank-3 solution of (18. ip which is implicitly defined by 

dh 

x, y) — C\(x — u x t) + — — — {x — u x t, y — u 2 t), 

or 1 

u 2 (t, x, y) = d(y - uH) ~^{x- uH, y - uH), ( 8 ' 13 ) 
a{t) = A 1 {{1 + dt) 2 + B.t 2 )- 1 ^, A, e R + , 

where the function h depends on two variables r 1 = x — w 1 t and r 2 = y — u 2 t and satisfies 
the nonhomogeneous Monge- Ampere equation 

h r i r ih r 2 r 2 — h r i r 2 = B\. (8.14) 

Depending on the selection of particular solutions of this equation we obtain Riemann double 
waves or other types of rank-2 solutions of (18.11) . 



9 Summary remarks 

The objective of this paper was to develop a new systematic way of constructing rank-k solu- 
tions of quasilinear hyperbolic systems of first order PDEs in many dimensions. Specifically, 
we have been interested in nonlinear superpositions of Riemann waves, which constitute the 
elementary solutions of these systems and are ubiquitous in the equations of mathematical 
physics. Interactions of Riemann waves are obviously present in many nonlinear physical 
phenomena. However there are still only a few examples of multiple rank solutions describ- 
ing them in multi-dimensional systems. Most of these solutions were obtained through the 
generalized method of characteristics. The main idea behind our approach has been to look 
at this type of solutions from a different point of view, namely, to reformulate them in terms 
of symmetry group theory. 

Let us now recapitulate our analysis. We look for rank-k solutions of the system (13.11) . 
expressible in terms of Riemann invariants, u = f(r l (x, u), . . . , r k (x, u)), where / : M k — > M. q . 
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Each Riemann invariant is associated with a specific wave vector involved in the interaction, 
i.e. r A (x, u) = Xf(u)x\ A = 1, . . . , k, where ker(Af A l {u)) ^ 0. The basic feature of these 
solutions is that they remain constant on (p — fc)-dimensional hyperplanes perpendicular 
to the set of linearly independent wave vectors A 1 , . . . , X k . In the context of group theory, 
this means that the graph {x,u(x)} of these solutions is invariant under all vector fields 
X a = £l(u)d x i with XfCa = for 1 < a < p — k. Then u(x) is the solution of (13. ip for some 
function /, because the set {r 1 , . . . , r k , u 1 , . . . ,u q } constitutes a complete set of invariants 
of the Abelian algebra L of such vector fields. The implicit form of these solutions leads 
to major difficulties in applying the classical symmetry reduction method to this case. To 
overcome these difficulties, we rectify the set of vector fields X a by a change of variables on 
X x U, choosing Riemann invariants as new independent variables. The initial equations (13.1 j) 
expressed in the new coordinates, complemented by the invariance conditions for the rectified 
vector fields X a , form an overdetermined quasilinear system (I3.2ip . Thus the solutions of 
this system are invariant under the Abelian group corresponding to L. The vector fields X a 
constitute the conditional symmetries of the initial system (13. ip . The consistency conditions 
for the overdetermined system (I3.2ip . that is, the necessary and sufficient conditions for the 
existence of conditionally invariant solutions of (13. ip . have been derived here and they take 
the form of the trace conditions (13.301) and (13.341) . Given these conditions, we were able to 
devise a specific procedure for constructing solutions in terms of Riemann invariants. We 
present it for the case of rank-2 solutions, however higher rank solutions can, in principle, 
be constructed by analogy. The computational difficulties should not be underestimated 
here and in many cases additional assumptions are needed in order to perform integration or 
to arrive at compact forms of solutions. Nevertheless, the implementation of the proposed 
CSM is still easier than of the GMC The latter imposes stronger restrictions on the wave 
vectors X A , which contribute to computational complexity as well as narrowing of the range 
of obtained solutions. 

As the application to the isentropic flow equations shows, our approach has proved quite 
productive. We were able to reconstruct the general rank-2 and rank-3 solutions obtained 
via the GMC and to deliver several new classes of solutions, namely in the cases E±Si, S1S2, 
EiE 2 Si, E1S1S2 and S1S2S3. For the equations of an isentropic flow with a sound velocity 
depending on time only (an assumption which simplifies things considerably) we obtained 
the arbitrary rank solution, together with the Cauchy conditions in a closed form. 

Moreover, we present a simple technique which allows us to overcome the main weakness 
of solutions expressible in terms of Riemann invariants, resulting from the fact that the first 
derivatives of Riemann invariants, in most cases, tend to infinity after some finite time. We 
show that a proper selection of the arbitrary functions appearing in the general solution 
can lead to bounded solutions even in the cases when Riemann invariants admit a gradient 
catastrophe. We obtained numerous such solutions which, to our knowledge, are all new 
(we include here only some of them, namely (RTTOj) . (l6"22l) . (jOT)) . fl633l) . fTT^j) . ([U]), dUSJ), 
(17.181) and (17.231) ). Most of these solutions have a soliton-like form and this fact is of note 
since the integrability properties of soliton theories do not easily generalize to more than two 
dimensions. 

Our technique is applicable to a very wide class of systems, which includes many physi- 
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cally meaningful models. Given the promising results obtained, we expect it may be useful in 
such areas as nonlinear field equations, general relativity or equations of continuous media. 
Let us note also that, though the notion of Riemann invariants was originally defined for 
hyperbolic systems only, it seems that it can be easily adapted to elliptic systems. Since 
the conditional symmetry method can be applied to these systems, it is worth investigating 
whether our approach to constructing rank-k solutions can be extended to the elliptic case. 
Some preliminary analysis suggests that to be feasible. 
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NSERC of Canada and FQRNT of Quebec. 
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